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Abstract. The QCD transition is studied on lattices up to N t — 16. The chiral 
condensate is presented as a function of the temperature, and the corresponding 
transition temperature is extracted. The equation of state is determined on lattices 
with Nt = 6, 8, 10 and at some temperature values with Nt = 12. The pressure and the 
trace anomaly are presented as functions of the temperature in the range 100 ...1000 
MeV . Using the same configurations we determine the continuum extrapolated phase 
diagram of QCD on the fx — T plane for small to moderate chemical potentials. Two 
transition lines are defined with two quantities, the chiral condensate and the strange 
quark number susceptibility. 
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Figure 1. Left: Subtracted chiral condensate A;. s as a function of the temperature. 
The gray band is the continuum result of our collaboration, obtained with the stout 
action. Right: The trace anomaly normalized by T 4 as a function of T on N t = 6, 8, 10 
and 12 lattices. The inset shows a comparison with the results of the Hadron Resonance 
Gas model, including resonances from the Particle Data Book up to 2.5 GeV mass. 



1. Introduction 

The study of QCD thermodynamics is receiving increasing interest in recent years. A 
systematic approach to determine the properties of the deconfinement phase transition 
is through lattice QCD. Lattice simulations indicate that the transition at vanishing 
chemical potential is merely an analytic crossover pQ. Some interesting quantities that 
can be extracted from lattice simulations are the transition temperature T c , the QCD 
equation of state and, for small chemical potentials, the phase diagram in the fi— T plane: 
we review the results on these observables that have been obtained by our collaboration 
using the staggered stout action with physical light and strange quark masses, thus 
m s/ m ud — 28 [21 For all details we refer the reader to Refs. [U El |6]. 

2. QCD transition temperature and Equation of State 

We present here the results for the chiral condensate, and extract the value of T c 
associated to this observable; for the values of T c obtained from other observables, which 
reflects the nature of the crossover transition, we refer the reader to Ref . [I] . The chiral 
condensate is defined as {ipip)q = T<91n Z/(dm q V) for q=u,d,s. It is an indicator for the 
remnant of the chiral transition, since it rapidly changes around T c . We calculate the 
quantity A J)5 , which is defined as [(0-0>«,t - W m s(^)<r]/[(^)z,o - rni/m s (xjjxjj) St0 ] 
for l=u,d. Since the results at different lattice spacings are essentially on top of each 
other, we connect them to lead the eye (see the left panel of Fig. [T|). The value of T c 
that we obtain from the inflection point of this observable is T c = 157(3) (3). 

Next we present our results regarding the equation of state; in the right panel of 
Figure [H the T dependence of the interaction measure is shown for the 2 + 1 flavor 
system. We have results at four different lattice spacings. Results show essentially 
no dependence on "a", they all lie on top of each other. Only the coarsest N t = 6 
lattice shows some deviation around ~ 300 MeV. On the same figure, we zoom in to 
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the transition region. Here we also show the results from the Hadron Resonance Gas 
model: a good agreement with the lattice results is found up to T ~ 140 MeV. 

In order to obtain the pressure, we determine its partial derivatives with respect to 
the bare lattice parameters, p is then rewritten as a multidimensional integral along a 
path in the space of bare parameters. To obtain the EoS for various m^, we simulate for 
a wide range of bare parameters on the plane of m u ^ and (3 (m s is fixed to its physical 
value). Having obtained this large set of data we generalize the integral method and 
include all possible integration paths into the analysis [HI [7] . We remove the additive 
divergence of p by subtracting the same observables measured on a lattice, with the 
same bare parameters but at a different T value. Here we use lattices with a large 
enough temporal extent, so it can be regarded as T = 0. 

3. The QCD phase diagram at nonzero quark density 

We provided results for the transition temperatures T c at vanishing chemical potential 
(/i = 0). Now we move out to the /i 7^0 plane (for the details see [6]). As one increases 
\i the transition temperature T c (/i 2 ) decreases. Let us parameterize the transition line 
in the vicinity of the vertical /i=0 axis as T c (fx 2 ) = T c (l — k ■ [i 2 /T 2 ). 

In order to determine the transition temperature as a function of [i we use two 
quantities which are monotonic in the transition region and do not depend on fi for 
zero or infinite temperatures. The transition temperature is defined as the temperature 
value at which these observables take their value as given by the inflection points of the 
curves at fi=0. 

The two observables we use are the renormalized chiral condensate and the 
normalized strange quark number susceptibility. In order to measure the /1 dependence 
of these quantities we apply reweighting. Since our lattices are quite large the full 
reweighting method is quite expensive. Therefore, we truncate the \x dependence of the 
weights at fi 2 order (this truncation of the original [U [9] method is usually called the 
Taylor method; for a recent application see [TO]). 

The strange quark number susceptibility is the second derivative of the partition 
function with respect to the strange chemical potential. It needs no renormalization, 
since it is related to a conserved current. It is useful to normalize it by T 2 , which 
provides a dimensionless combination. It is easy to see that at T = its value is 0, 
whereas for infinitely large temperatures it approaches the Stefan-Boltzmann limit of 1. 

For the chiral condensate we apply here a slightly different renormalization 
prescription than what was used for the determination of T c . We cancel the additive 
divergences by subtracting the T = contribution, while the multiplicative divergence 
due to the derivative with respect to the mass can be eliminated with a multiplication 
by the bare quark mass. Then, in order to have a dimensionless combination the whole 
expression can be divided by the fourth power of some dimensionful mass scale. In 
this work we use the T=0 pion mass for the normalization. This observable has also \x 
independent limiting values at zero and at infinitely high temperatures (at T=0 this is 
true for chemical potetials smaller than the baryon mass, well within our applicability 
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Figure 2. Left: the pressure normalized by T 4 as a function of T on N t = 6, 8 and 10 
lattices. Right: The crossover transition between the cold and hot phases is represented 
by the coloured area (blue and red correspond to the transition regions obtained from 
the chiral condensate and the strange susceptiblity, respectively) . The lower solid band 
shows the result for T c (/i) defined through the chiral condensate and the upper one 
through the strange susceptibility. The width of the bands represent the statistical 
uncertainty of T c (n) for the given /x coming from the error of the curvature for both 
observables. The dashed line is the freeze-out curve from heavy ion experiments. The 
center of mass energies of these experiments are also shown. 



region) . 

Our final result is shown in the right panel of Figure EJ The crossover regions 
extent changes little as the chemical potential increases, and within it two definitions 
give different curves for T c (/i). It is useful to compare the whole picture to the freeze-out 
curve which summarizes experimental results on the T-/x points where hadronization of 
the quark-gluon plasma was observed. This curve is expected to lie in the interior of 
the crossover region, as is indicated by our results as well. 
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